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Summary . 

The  problem  of  a  semi-infinite  mode  III  crack  that 
suddenly  begins  to  propagate  at  a  constant  speed  Is 
considered  for  a  general  linear  viscoelastic  body.  A  simple 
closed  form  expression  for  the  Laplace  transform  of  the 
energy  release  rate  (ERR)  is  derived  under  the  assumption 
that  a  Barenbatt  type  failure  zone  exists  at  the  crack  tip. 
The  first  two  terms  of  a  short  time  asymptotic  series  for 
the  ERR  is  constructed  and  the  rate  at  which  the  ERR 
converges  to  steady-state  is  studied.  It  is  shown  that  the 
rate  of  convergence  to  steady-state  is  dependent  upon  crack 
speed  and  material  properties.  Moreover,  it  is  found  that 
whether  or  not  a  failure  zone  is  incorporated  into  the  model 
significantly  influences  both  quantitatively  and 
qualitatively  the  short  and  long  time  behavior  of  the  ERR. 
This  difference  is  important  to  predictions  of  stable  vs 
unstable  crack  speeds  based  upon  a  critical  ERR  fracture 


criterion. 
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1.  Introduction. 

Several  analytical  studies  of  dynamically  propagating 
cracks  in  linear  visoelastic  material  have  appeared  in  the 
literature  since  Willis (1967)  presented  an  analysis  of  the 
dynamic,  steady-state  propagation  of  a  semi-infinite,  mode 
III  (anti-plane  shear)  crack  in  an  infinite  viscoelastic 
body.  Employing  transform  methods  and  the  Wiener-Hopf 
technique,  Willis  constructed  the  dynamic  stress  intensity 
factor  (SIF)  for  a  special  class  of  crack  face  loadings  and 
for  a  standard  linear  solid  material  model.  Subsequently, 
Atkinson  and  List (1972)  introduced  transient  effects  into 
the  problem  by  assuming  that  the  crack,  initially  at  rest, 
begins  to  propagate  at  a  constant  speed  under  the  action  of 
suddenly  applied  loads  on  the  crack  faces.  Also  utilizing 
the  Wiener-Hopf  method,  they  derived  an  expression  for  the 
Laplace  transform  of  the  time  dependent  SIF  from  which  long 
and  short  time  asymptotic  approximations  and  numerical 
Laplace  inversion  calculations  were  obtained.  However, 
their  analysis  was  limited  to  consideration  of  constant 
applied  load  along  the  crack  faces  and  the  required 
Wiener-Hopf  factorization  was  effected  only  for  a  Maxwell 
material  model  and  the  Achenbach-Chao( 1962 )  three  parameter 
approximation  to  the  standard  linear  solid. 


Somewhat  later,  Atkinson  and  Coleman( 1977)  used  a 
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matched  asymptotic  expansion  technique  to  develop  an 
approximate  analysis  of  the  steady-state  propagation  of  a 
semi-infinite  mode  I  (plane  strain  or  plane  stress)  crack: 
propagating  in  a  clamped  viscoelastic  strip.  Shortly 
thereafter,  Atkinson( 1979 )  presented  an  approximate  analysis 
of  the  mode  I  counterpart  to  the  mode  III  problem  considered 
by  Atkinson  and  List(1972).  Their  argument,  ostensibly 
valid  for  fairly  general  material  models,  involved  slightly 
modifying  the  exact  elastic  result  of  Baker(1962)  in  order 
to  approximate  the  Laplace  transform  of  the  actual  dynamic 
viscoelastic  SIF.  The  dominant  term  for  each  of  the  short 
and  long  time  asymptotic  expansions  of  this  approximate  SIF 
was  derived  for  each  of  three  different  applied  crack  face 
loads:  a  constant,  a  delta  function  point  force  and  an 

exponentially  decreasing  form.  Also  in  that  paper,  Atkinson 
reconsidered  the  mode  III  problem  and  extended  the  Atkinson 
and  List  analysis  to  handle  the  above  three  types  of  crack 
face  loadings.  However,  consideration  was  limited  again  to 
the  Achenbach-Chao  material  model.  Atkinson  also 
constructed  an  expression  for  the  energy  release  rate  (ERR) 
based  upon  a  local  (i.e.  at  the  crack  tip)  work  argument  and 
the  singular  stress  field. 

Also  in  that  same  year,  Atkinson  and  Popelar(1979) 
presented  an  analysis  of  the  transient  constant  crack  speed 
mode  III  problem  for  a  viscoelastic  strip.  Constitutive 


relations  in  terms  of  differential  operators  were  assumed 
and  the  external  load  consisted  either  of  constant 
displacement  of  the  upper  and  lower  layer  boundries  or 
constant  tractions  on  them.  The  crack  faces  were  assumed  to 
be  stress  free.  Again  the  Wiener-Hopf  method  was  used  to 
construct  an  exact  expression  for  the  Laplace  transform  of 
the  SIF.  The  required  Wiener-Hopf  factorization  was  carried 
out  modulo  a  term  involving  a  Cauchy  type  integral . 

Atkinson  and  Popelar  then  restricted  attention  to 
numerically  approximating  the  Cauchy  integral  for  the 
steady-state  limit  case  and  assumed  a  standard  linear  solid 
material  model . 

A  year  later,  Atkinson  and  Popelar ( 1980)  addressed  the 
more  difficult  mode  I  problem  for  a  viscoelastic  strip. 

Again  by  use  of  the  Wiener-Hopf  method,  a  formal  expression 
for  the  Laplace  transform  of  the  SIF  was  constructed 
containing  a  complicated  Cauchy  integral.  As  with  the 
corresponding  mode  III  case,  the  integral  was  studied 
numerically  in  the  limiting  special  case  of  steady-state 
crack  propagation  in  a  standard  linear  solid. 

Somewhat  later,  Walton(1982)  examined  further  the 
steady-state  mode  III  problem  considered  by  Willis( 1967) . 
Utilizing  the  Riemann-Hilbert  rather  than  the  Wiener-Hopf 
methodology,  he  constructed  a  simple  closed  form  expression 
for  the  SIF  valid  for  general  crack  face  loadings  and  very 
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general  material  models.  More  specifically,  constitutive 
equations  expressed  in  terms  of  convolution  integrals  rather 
than  differential  operators  were  adopted  and  the  results 
were  shown  to  be  valid  irrespective  of  any  assumed  time  rate 
of  decay  of  the  viscoelastic  shear  modulus.  In  contrast, 
constitutive  relations  in  terms  of  constant  coefficient 
differential  operators,  necessarily  force  an  exponentially 
decaying  modulus  thereby  preventing  consideration  of  the 
important  class  of  power-law  models  which  more  effectively 
represent  the  mechanical  response  of  many  real  viscoelastic 
materials,  such  as  rubber,  than  do  exponentially  decaying 
functions . 

Subsequently,  Walton(1983)  extended  the  above  analysis 
to  determine  the  angular  dependence  of  the  stress  field  in  a 
neighborhood  of  the  crack  tip.  In  particular,  it  was  shown 
that  the  asymptotic  stress  field  at  the  crack  tip  has  the 
same  angular  dependence  as  the  corresponding  dynamic  elastic 
problem.  Only  the  SIF  differs  between  the  elastic  and 
viscoelastic  fields. 

Walton(1985)  next  considered  the  steady-state  mode  III 
problem  for  a  viscoelastic  strip.  Again  utilizing  the 
Rlemann-Hilbert  method,  a  closed  form  expression  for  the  SIF 
was  constructed  for  general  loadings  and  shear  modulus.  The 
form  of  the  solution  exhibits  clearly  the  combined  effects 
of  material  properties,  crack  speed  and  layer  thickness  upon 
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the  SIF. 

More  recently,  Walton( 1987a)  reconsidered  steady-state 
mode  III  crack  propagation  In  an  Infinite  viscoelastic  body 
in  order  to  investigate  the  implications  of  including  a 
failure  zone  model  of  Barenblatt ( 1962 )  type  into  the 
determination  of  the  ERR  from  a  global  energy  balance 
calculation.  A  simple  closed  form  expression  for  the  ERR, 
which  in  this  case  is  just  the  work  done  by  the  tractions  in 
the  failure  zone,  was  derived  under  the  same  mild  conditions 
on  the  shear  modulus  assumed  in  Walton(1982)  and  for  a 
fairly  broad  class  of  crack  face  and  failure  zone  loadings. 
It  was  then  observed  that  whether  or  not  a  failure  zone  is 
incorporated  into  the  model  greatly  influences  both 
qualitatively  and  quantitatively  the  dependence  of  the  ERR 
upon  crack  speed  and  material  properties.  In  particular, 
calculations  based  upon  a  failure  zone  seem  to  reflect  more 
closely  experimental  observations  of  cracks  rapidly 
propagating  in  real  viscoelastic  material.1 

The  methods  of  Wal ton ( 1987a)  can  be  applied  to  the 
calculation  of  the  ERR  for  a  wide  variety  of  dynamic 
viscoelastic  crack  problems.  Schovanec  and  Walton( 1987c ) 
recently  completed  the  analysis  of  the  dynamic  steady-state 
propagation  of  two  parallel  mode  III  cracks  in  an  infinite 

1  Private  communication  with  Prof.  W.G.  Knauss  and  Prof.  K. 


Kuo . 
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viscoelastic  body.  Also,  Walton{ 1987b)  has  recently 
completed  a  study  of  the  mode  I  analog  of  Walton( 1987a ) .  In 
both  of  these  investigations  a  Barenblatt  failure  zone  model 
was  adopted.  Of  related  interest  are  two  additional  papers 
by  Schovanec  and  Walton( 1987a, 1987b)  in  which  these  same 
methods  were  applied  to  quasi-static  mode  I  crack 
propagation  in  non-homogeneous  viscoelastic  material.  It 
should  also  be  noted  that  Knauss(1973)  and  Schapery ( 1975 ) 
applied  the  Barenblatt  model  to  quasi-static  viscoelastic 
crack  growth  in  viscoelastic  material  and  observed  that 
whether  or  not  a  failure  zone  is  incorporated  into  the  model 
greatly  affects  the  behavior  of  the  ERR. 

The  present  paper  applies  the  above  Barenblatt  failure 
zone/Riemann-Hilbert  program  to  the  transient  mode  III 
problem  considered  by  Atkinson  and  List (1972).  In 
particular  a  simple  closed  form  expression  for  the  Laplace 
transform  of  the  ERR  is  derived  for  a  very  general  class  of 
material  models.  Only  crack  face  loadings  of  a  simple 
exponentially  decaying  type  are  considered  but  the  extension 
to  more  general  loadings,  described  in  Walton ( 1987a) ,  offers 
no  difficulty  in  the  present  context  and  is  omitted  for  the 
sake  of  brevity.  Moreover,  the  loadings  assumed  here 
contain  all  of  the  essential  ingredients  of  the  Barenblatt 
model  and  as  such  are  amply  suited  to  Illustrating  the 
qualitative  features  expected  from  more  general  mathematical 
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models . 

From  the  Laplace  transform  of  the  ERR  both  long  time 
and  short  time  asymptotic  analyses  are  considered.  In 
particular,  the  first  two  terms  of  a  short  time  asymptotic 
series  is  derived  and  the  question  of  the  convergence  rate 
to  the  steady-state  solution  is  addressed.  It  is  shown  that 
the  rate  of  convergence  to  the  steady-state  solution  is 
greatly  dependent  upon  crack  speed  and  material  properties. 
Moreover,  whether  or  not  a  failure  zone  is  assumed  to  exist 
greatly  influences  both  the  short  time  and  long  time 
behavior  of  the  ERR  as  well  as  predictions  of  stable  vs 
unstable  crack  speeds  based  upon  a  critical  ERR  fracture 


criterion. 
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2.  Problem  Formulation  and  Stress  Analysis. 

The  problem  to  be  considered  is  that  of  a  semi-infinite 
mode  111  crack  that  begins  to  propagate  at  a  constant  speed 
v  in  an  infinite  viscoelastic  body  due  to  the  sudden 
application  of  crack  face  tractions  that  then  travel  with 
the  crack.  The  governing  field  equations  for  the  motion  of 
a  linear  viscoelastic  solid  are 

where  o^,  and  ui  denote  the  stress,  strain,  and 

displacement  fields  respectively.  In  (2.1),  p* dt  denotes 

the  Riemann-Stielt jes  convolution  p*dt  =  p(t-r)  de(r). 

J  —00 

Since  the  deformation  is  assumed  to  be  antiplane 
strain,  ul=0,  ua=0,  and  the  only  equation  of  motion  not 
identically  satisfied  is  p*cLiu1=pu,  . 

A  semi-infinite  crack  lying  along  the  negative  Xj-axis 
is  assumed  to  begin  to  propagate  at  time  t=0  with  a  constant 
speed  v  driven  by  loads  o t , (x, , 0 , t ) =f (x, -vt )  which  follow 
it.  The  corresponding  initial-boundary  value  problem  is 


pu,=p*ddu,  =A 

p(0)u,  (x,  ,: 

t,+j 

u,  (X,  ,X2  ,T  )  p' 

0 

(t- 

r)  dr 

(2. 

.2) 

with  initial 

conditions 

u. 

=0, 

u,*0  at  t=0 

(2. 

.3) 

and  boundary 

conditions 

e2 

t  (Xj 

,0,t)®f (Xj-vt) 

Xx 

<  vt 

(2. 

.4) 

u 

»  (X, 

,0,t)=0 

X, 

>vt 

»  (Xt 

2 

,x2  ,  t  )-*0  as  xA 

+x, 

2^. 
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From  (2.1)  it  can  be  seen  that 
On  (x,y,t)=|^  (p* du,  ) 

a  r  r* 

p(0)u,  (x,  ,x2  ,t)  +  u,  (x,  ,x2  ,T  )p*  (t-r  )  dr  . 

*  L  Jo 

It  is  convenient  to  change  from  the  fixed  coordinates 
(x, ,x2 ,t)  to  the  moving  coordinate  system  (x,y,t)  given  by 
x=Xi~vt,  y=x2  and  to  define  w(x,y,t)  by 


w(x,y , t ) =w(Xj -vt , x2 , t ) =u, (xt , x2 , t ) .  In  the  moving 
coordinates,  equation  (2.2)  becomes 

rd  a  1 2  f  r* 

p  3T"VJx  w(x»Y»t)ss<4|>(0)w(x,y.t)+j  w(xt  -vr  ,x2  ,r  )p*  (t-r  )  dr 

-A  p (0)w(x,y, t )+ [  w(x+vt-vr  ,y,r  )p‘  (t-r  )  dr  . 

Jo 

Application  of  the  Fourier  transform,  defined  by 
f(p,Y,t)=J  eipxf(x,y,t)  dx,  to  equation  (2.6)  results  in 


ra 

p[at+ivp] 


w(p,y,t)= 


|^-ip)2+— 2j  jp(0)w(p,y,t)+J  w(p,y,r  )e”iVp(  1  r)p’(t-r)  dr 
Subsequent  application  of  the  Laplace  transform 


1 » Y » s )  =  J 


9(P»Y>r)e  dr  to  (2.7)  yields  the  equation 


p(s+ivp) 2w(p,y,s)=^— 2-p2| ^p (s+ivp)w(p,y,s) 

in  which  p(s)  is  the  Carson  transform  of  thfe  shear  modulus 


given  by  p(s)=sp(s) 


=A*  ( 0 )  +  J 


-sr  . 
e  dp . 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


Equation  (2.8)  can  be  rewritten  as 
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2  -  _ 

— _w-  [ p2+- — ^ - (s+ivp)2lw=0  whl 

ay^  *•  P(s+ivp)  J 

w(p,y,s)=A(p,s)e~^s'p)  lyl 

where  fl(s,p)=  p2+- — - - (s+ivp)2| 

*•  /S(s+ivp)  ■* 


which  has  the  solution 


(1/2) 


must  be  chosen  so 


that  Re  p>0. 


In  a  similar  manner,  the  Fourier  and  Laplace  transforms 
may  be  applied  to  the  constitutive  equation  (2.5)  to  produce 

On  (P»y.s)=p(s+ivp)|^w(p,y,s)  . 

oo 

If  one  defines  f+(p)  and  f  (p)  by  f+(p)=[  eipxf(x)  dx  and 

J0 

.0 

f  (p)  =  |  elpxf(x)  dx,  then  the  boundary  condition  (2.4) 

J  —CO 

transforms  to 

*  ^  ■“  _  ^  ^  * 
o,»(p.0,s)+oa,  (p,0,s)=f(  (s+ivp)y^w(p,  0,s) 

=-/3{s,p)?(s+ivp)w(p,0,s)  . 

From  (2.4)  it  can  be  seen  that  o2t(p,0,s)=f(p,s)  and 
w(p,0,s)=w  (p,0,s).  It  is  assumed  a  priori  (and  is  easily 


(2.9) 


(2.10) 


(2.11) 


verified  a  posteriori)  that  w  (p,0,s)  and  oaa( p,0,s)  have 

analytic  extensions  w-(z,0,s)  and  oa*( z,0,s)  for  lm(z)<0  and 

lm(z)>0,  respectively,  which  vanish  as  |z (-*«►.  Thus  the 

transformed  boundary  condition  (2.11)  can  be  recast  as  the 

Riemann-Hilbert  problem:  find  F+(z)  analytic  for  Im  <z)>0 

and  F  (z)  analytic  for  lm(z)<0  such  that 

lim  F+(z)»  lim  F  (z)=0  and  on  Im(z)B0, 

Ia(z)-»+»  Im(z)-»-«> 
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F+(p)=T(p)F-(p)-f (p,s)  for  pe  ( -»,») 


where  F+(z)=ott(z,0,s),  F~(z)=w  ( z ,  0 ,  s )  ,  and 

r  2  p  2i(1/2) 

T(p) =-p ( s+ivp )  p^+- — - - (s+ivp) 

1  fJ  ( 8  +  1  Vp )  ■> 


It  is  well  known  that  the  solution  of  (2.12)  is 


F±(z)=X±(z)^iIJ 

*  aj-~  X  (r)  r  s 


where  X  (z)  solves  the  homogeneous  Riemann-Hilbert  problem 
X+(p)=T(p)X-(p)  for  p*.  (-<*,<*>)  . 

To  solve  (2.14),  it  is  convenient  to  factor  T(p)  into 


the  product 


T ( p ) -Gt (p)Ga ( p ) G, (p) 

r  is  91(1/2) 

in  which  G, (p)=-p(s+ivp) ,  Ga (p)*[(p--^)  J  ,  and 

2  2  ( 1  /  2 ) 

G, (p)=[ — He  b  ~  — -  1  •  X*(z)  may  now  be 

l(p-=5)  P (s+ivp) J 

constructed  as  the  product  X* (z)=X,± (z)X2± (z)X,* (z)  with 
each  X1± (z)  satisfying  the  Riemann-Hilbert  problem 

Xi+(p)=Gi(p)Xi”(p) • 

What  will  ultimately  be  required  is  to  solve  the 
homogeneous  Riemann-Hilbert  problem  (2.14)  for  each  fixed  s 
on  a  Bromwich  path  with  Re(s)>0.  This  will  be  accomplished 
by  first  assuming  8  real  and  positive  and  then  invoking  an 
analytic  continuation  argument.  Additionally,  for  the 
subsequent  analysis  the  shear  modulus  will  be  assumed  to  be 
positive,  continuously  differentiable,  non-increasing. 


(2.12) 


(2.13) 


(2.14) 


(2.15) 
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convex  and  such  that  />(«*)  *lim  fj(  t)>0.  Convexity  is 

t-*«o 

sufficient  but  certainly  not  necessary  to  insure  the 
validity  of  the  following  calculations  and  though 
theoretically  overly  restrictive,  it  holds  for  most  of  the 
customary  models  such  as  a  standard  linear  solid  or  a 
power-law  material.  Moreover,  it  is  worth  noting  that  no 
explicit  time  decay  rate  for  the  shear  modulus  needs  to  be 
specified  for  the  results  to  be  valid.  From  the  fact  that 
fj(  t)=0  for  t<0,  it  easily  follows  that  [p  ( s+ivz)]-1  is 
analytic  for  Im(z)<0.  Therefore  one  may  choose 

X, +(z)=l  and  Xj'tz)*-!? (s+ivp) ]-1 .  (2.16) 

Since  the  product  Ga  (p)Ga  (p)*/J(s,p)  ,  the  branches  of  Ga(p) 

and  Ga(p)  must  be  chosen  so  that  their  product  satisfies  the 

requirement  that  Re  /!(s,p)>0.  This  condition  can  be  met  by 

choosing  the  branch  of  z^1/<2^  with  branch  cut  along  the 

negative  real  axis  for  both  Ga(p)  and  G,(p).  (See 

figure  2.1.)  Therefore  Ga  (p)  can  be  expressed  as 
is  1 

Ga  (p)=sgn(p)  (p-— )  and  Xa  (z)  may  be  chosen  to  be 

Xa+(z)=»+(z)  and  Xa"(z)*w"(z) (p-l|)_1  (2.17) 

in  which  w+(z)*z^ 1/'2)  with  branch  cut  along  the  negative 
Imaginary  axis  and  w~  (z)=z^ with  branch  cut  along  the 
positive  imaginary  axis.  Finally,  one  may  construct  Xa1(z) 
by 

Xa±(z)=exp (r±(z))  where  '■‘(zJ-jiLf  l0?^»(r))  dr. 

"  —00 


(2.18) 
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To  find  a  closed  form  expression  for  Xt±(z),  it  is 
necessary  to  determine  the  mapping  properties  of  G»(z)  for  z 
in  the  half-plane  lm(z)<0.  If  one  first  considers 
p(s+ivz)*p( s-q+ivp)  on  the  horizontal  lines  z=p+iq,  q<0, 

(-«,«•),  then  it  follows  easily  from  the  stated  assumptions 
on  p(  t)  that 

(i)  P ( 0)  =p  («*)  <p ( s-q)  <Re  p  (s-q+ivp)  <p  ( 0)=p  (s-q±  i<»)  ; 

(ii)  Im  p ( s-q+ivp ) =-Im  p(s-q-ivp); 

(iii)  arg  p( s-q+ivp) ; 

(iv)  lim  p (s-q+ivp )=p ( 0 )  with  p(s-q)  converging 
monotonically  to  p(0) . 

Therefore  p(s+ivz)  maps  these  lines  to  the  curves  shown  in 
fig  2.2a. 

The  linear  fractional  transformation  S ( z)  =z_j^8/>V)  maps 
the  lines  z=p+iq,  q<0,  pM-®9,®9)  to  the  circles  shown  in 
fig  2.2b  and  S(p+iq)  exhibits  the  following  properties: 

(v)  0=S(0)<S(0+iq)<Re  S(p+iq)<S(±«o+iq)=l : 

(vi)  Im  S(p+iq)=-Im  S(-p+iq); 

(vii)  arg  S(p+lq)|j°  ; 

(vlii)  lim  S(p+iq)*l  with  S(iq)  converging 

monotonically  to  1. 

If  the  elastic  shear  wave  speeds  corresponding  to 

2 

infinite  and  zero  time  are  defined  by  c„  «^(«»)/p  and 
c2«p(0)/p  then  from  (i)-(viii),  it  follows  that  [G,(z)]2  has 


16 


the  following  properties: 

(ix)  Im  G, 2{p+iq)=-Im  G,2(-p+iq); 

(x)  arg  G,2(p+iq)j^  ; 

(xi)  lim  G, 2 (p+iq) =1- ( v/c ) 2 ; 
q-»-«o 

(xii)  -(v/c,)2=G,2(0)<G,2(iq)<G,2(-i<»)  =  l-(v/c)2  where 

2  2 
G,  (iq)  is  monotonically  increasing  to  l-(v/c) 

as  q-»-«> . 

2 

Thus  it  can  be  seen  that  G,  (z)  maps  the  horizontal  lines 

z=p+iq,  q<0,  pe  (-<»,«>)  to  the  curves  shown  in  fig  2.2c. 

2 

Furthermore,  it  can  be  seen  from  (xii)  that  G,(z)  has  a 
unique  root  z*=iq,  in  the  half-plane  lm(z)<0  for  any 
positive  real  value  of  s. 

If  the  branch  cut  for  the  square  root  defining  G,(z)  is 

chosen  along  the  negative  real  axis  then  G,(z)  is  analytic 

for  lm(z)<0  except  for  the  branch  cut  on  a  segment  of  the 

negative  imaginary  axis  across  which  it  has  a  jump 

discontinuity  given  by  G,  (±0+iq)=±l  | (syv^.qj  for  <3*<(3^°* 

It  follows  that  log  G,(z)  is  analytic  for  lm(z)<0  away  from 

this  line  segment  and  has  the  jump  discontinuity 

log(G, (+0+iq) )-log(G, (-O+iq) )»in  across  z=iq,  q*<q<0.  One 

00 

can  then  evaluate  r±(z)»jA^|  2 2  <jr  as  in 

Walton(1982)  and  conclude  that 

r* (z)«( 1/2) ^log(q,i-z)-log(-z)+log(G, (z) ) 

fO  for  lm(z)>0  ] 

\-log(G,(z))  for  lm(z)<0  J' 
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Therefore  (2.18)  reduces  to 
X,+(z)»iw+(z-iq,)[G,(«)](1/2)/u+(z) 

X,“(z)=iw-(z-iq,)IG, («)](1/2)/w'(z)G, (z). 

From  equations  ( 2 . 16) - ( 2 . 19 )  it  then  follows  that 
X+(z)=iw+(z-iq,)(G, («)](1/2) 

X*(z)  =-!»'{ z-iq,)[G,(«*)](  1/2  )/[?{s+ivz)  (z-i(s/v)  )G,  (z)  )  . 
Finally,  for  a  specific  load  o~(x)=f(x)  one  can 

determine  F+(z)=o2*(z,0,s) ,  F_ ( z ) =w- ( z , 0 , s )  from  equations 
(2.13)  and  (2.20).  The  Laplace  transform  of  the  stress 
intensity  factor  (SIF)  K(s)  can  now  be  calculated  as  in 
Walton( 1982) .  In  particular,  it  is  straight  forward  to  show 
that 


o(x)~5!^.x(-1/2)  as  x-*0+ 

-w _ /^,(l/2)e(~1/4)”1r  f{r 


where  K(s)  =  lG,  («»)  ] 


J  lill 

J-oo  x+(r) 


Again,  as  in  Walton( 1982 ) ,  (2.21)  may  be  simplified  to 


K(s)^iJ~  li£)|*|(‘1/2,e"Xq*  dx. 


(2.19) 


(2.20) 


(2.21) 


(2.22) 


in  which  it  should  be  emphasized  that  q,  is  a  function  of  s 
determined  implicitly  by  G»(iq„)*0. 
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3.  Calculation  of  the  Energy  Release  Rate. 

The  Energy  Release  Rate  (ERR)  will  now  be  calculated 

based  upon  the  assumption  that  a  Barenblatt  type  failure 

zone  exists  at  the  crack-tip.  Specifically,  it  is  assumed 

that  two  loads  are  acting  on  the  crack-faces:  the  applied 

(external)  tractions  denoted  og(x)  and  the  cohesive 

(failure)  stresses  o~(x)  acting  in  a  failure  zone  of  length 

af  immediately  behind  the  crack-tip.  The  only  assumptions 

about  o~(x)  are  that  a^  is  small  relative  to  some  length 

scale  a  associated  with  o”(x)  and  that  K  +K  =0  where  K  and 
e  e  e  f  e 

Kf  are  the  SIF's  corresponding  to  a~  and  a~,  respectively. 
Hence  the  effect  of  the  failure  zone  is  to  cancel  the 
singular  stresses  ahead  of  the  crack-tip  and  thereby  produce 
a  cusp  shaped  crack  profile  behind  the  tip. 

The  ERR,  G(t)  (defined  to  be  the  energy  flux  into  the 
crack  tip  per  unit  crack  advance)  is  given  by 
1  rvt 

G(t)«-  of (xj-vt)u, (x,-vt,0,t)dxj 

J vt-af 

which  in  the  moving  coordinates  becomes 

G(t,4j  [lt‘v^]w(x'°'t)dx- 

"af 

It  is  impractical  to  calculate  (3.1)  directly  since  the 
expression  £^-v|^jw(x,0,t)  has  a  verY  complicated 
dependence  on  o@(x)  and  of(x)  through  the  inverse  Laplace 


(3.1) 
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transform  of  (2.13)  and  (2.20).  A  computationally 
convenient  expression  will  be  derived  for  the  Laplace 
transform  of  G(t),  G(s),  for  the  special  loadings  a  and  o 

0  I 

used  in  Walton( 1987a) .  Specifically,  it  is  assumed  that 

o“(x)=Leexp(x/ae)H(-x)  and  a “ (x) =-Lfexp(x/af )H( -x)  (3.2) 

where  a. /a  <<1.  For  a. /a  small  enough,  the  fact  that  o~(x) 
re  x  e  r 

does  not  have  support  in  some  small,  compact  interval  behind 
the  crack-tip  will  have  a  negligible  effect  on  the  results. 

The  assumptions  (3.2)  clearly  incorporate  the  salient 
features  of  the  Barenblatt  model,  namely  an  applied  load  o~ 
paired  with  a  cohesive  stress  o~,  each  with  an  associated 
length  scale  afi  and  af  respectively,  such  that  a~  cancels 
the  singular  stress  produced  by  a~  and  af/ae<<l.  It  should 
be  noted  that  in  this  case  (3.1)  should  be  replaced  by 

o  r 

G(t)=H  °f(x)  w(x,0,t)dx.  (3.3) 

^  ■—OO  l 

As  shown  in  Walton  (1987a),  it  is  straight  forward  to  extend 
the  analysis  to  treat  more  general  loads  in  the  form 

CO 

a  (x)*L|  etx>/adh(t)  where  h(t)  is  any  signed  measure  for 

Jo 

which  the  integral  makes  sense.  However,  for  the  sake  of 
brevity  that  development  is  not  Included  here. 

If  one  applies  Parseval's  formula  and  then  takes  the 
Laplace  transform  of  the  resulting  expression,  the  Laplace 
transform  of  the  ERR  is  found  to  be  given  by 
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G(s)=iJ  o~(p) |s+ipvjw"(p,0,s)dp 


(3.4) 


where  of{ p)  denotes  the  inverse  Fourier  transform  of  o~(x). 

For  the  particular  choice  of  o“(x)  given  in  (3.2), 

v_  v_ 

o f  ( P ) ja^  <j }  .  It  is  easily  seen  that  o"(p)  has  a 

v_ 

meromorphic  extension  o~( z)  with  a  simple  pole  at  z=-i/af. 

Furthermore,  since  w~( z , 0 , s ) =F-( z )  is  analytic  for  Im(z)<0, 
(3.4)  can  be  evaluated  by  residues  and  whence  it  follows 
that 

"kf 

G(s)=-^-(s+(v/af )  )F  (-i/af). 

It  remains  to  evaluate  F  (~i/af).  To  this  end,  one  begins 
by  first  noting  that  (2.11)  can  be  rewritten  as 


(3.5) 


F  (p)=o (p)/T(p) 


(3.6) 


where  o(p)*o  (p)+o  (p). 

From  the  Barenblatt  model  o”(x)=o’(x)+oj(x)  and 
consequently 


(3.7) 


o~(P)-°~(P)+of (P) • 


(3.8) 


Also  o  (p)  in  (3.7)  can  be  determined  by  application  of  the 
Pleaelj  formula  to  (2.13)  thereby  obtaining 


i  +  (P)»F+(p)— Jf(p)+^iTX+(p)[  -2l. 

2  2nl  X+(r)  P 


This  can  be  rewritten  using 


(3.9) 


o  +  (P)-Og(p)+flJ(p) 


(3.10) 
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wherein  o^(p)  and  o*(p)  are  given  by 


J -»  X  ( r  )  e 


From  (3.2)  it  can  be  seen  that  °e(P)=is(p-( ±/a  ) )  and 


(3.11) 


therefore  a  ^(z)  is  meromorphic  with  pole  z=i/ag.  Since 
X+(z)  is  analytic  for  lm(z)>0,  (3.11)  can  be  calculated  by 
residues  with  the  result 

«j(p)~ «;<p>[i  -  V*1 


X' (i/ae) 


Similarly, 


’ ^  ( P)  —~o  j  ( p )  f  1  -  *  lV]-  | 
r  r  L  X+ ( i/a, ) J 


(i/af) 

which  combined  with  (3.10)  and  (3.12)  yields 


o(P)-^[ - 

as  l(p-(i 


(p-(i/ae) )X+(i/ae)  (p-(i/af ) )X+(i/af) 
Equation  (3.14)  can  now  be  simplified  by  the  Barenblatt 
hypothesis  Ke+Kf=0.  For  the  special  loadings  given  by 
(3.2),  equation  (2.21)  yields 


and 


Kf(s,.CG,C)]{1/2)e<-1/^"i 


sX  (i/ae) 

-L. 


,  which  under  the 


sX  ( i/af ) 

» 

Barenblatt  hypothesis  results  in  the  identity 


(3.12) 


(3.13) 


(3.14) 


X+(i/ae)  X+(i/af) 

Equations  (3.14),  (3.15),  (3.6)  and  the  observation 


(3.15) 
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that  T(p)=X  (p)/X~(p)  then  yield  that 
a  L  (a  -a  )  x"(-i/a.) 

F  (-1/af>=£-Tsf  Tir+iT)  — - 1  • 

zs  iae+af)  x  (i/ae) 

One  may  now  substitute  (3.15)  and  (3.16)  into  (3.5)  thereby 


obtaining 

L  2  a. (a  -a  )  X-(-i/a  )X+(i/a  ) 

G(s)^is+{v/af  1)-^) 


(X+(i/ae)]2 


which  when  combined  with  the  result  (2.20)  derived 

previously  for  X±(z)  produces  the  desired  expression  for  the 

Laplace  transform  of  the  ERR 
2. 


L  ‘a  (a  -af) 

G(s)s?5  ~ra-^) 


2, (1/2) 


pv 


(l+saf/v)  p(s+v/af) 


~(l/2) 


[l-aeq,]  p(s+v/af) 

In  order  to  present  the  results  in  a  nondimensional 

form,  it  is  necessary  to  introduce  certain  parameters. 

First,  a  nondimensional  shear  modulus  is  defined  by 

p(t)=p  m(t/r)  where  p  =lim  p(t)  and  thus  lim  m(t)=l.  Also, 

t~»°°  t-» 


the  nondimensional  parameters  y,  t,  a,  and  p  are  defined  by 
”»!Bv/c*I  *«af/ae,  o=c,r/ae,  and  /3=-aeq„  where  0<-»<c/c„,  e<<l, 
o>0,  and  p>0. 

The  Laplace  transform  of  the  ERR  can  now  be  rewritten 
in  nondimensional  form  as 


(3.16) 


(3.17) 


(3.18) 
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sG(s)~ 


Le2% 


U-*) 

TTf*) 


U-(M)2Iu/2) 
TI+?)  ' 


(a/fe  )2 
(rs+ai/e  )  2 


1 

m(r s+aK /t ) 


1 


m( r s+a? /«. ) 
-U/2) 


It  is  the  study  of  this  basic  expression  for  G(s)  found  in 
(3.19)  that  leads  to  the  special  cases  and  asymptotic 
approximations  found  in  the  next  section. 


(3.19 
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4.  Short-time  and  Long-time  Asymptotic  Solutions. 

Asymptotic  expansions  for  the  ERR,  G(t),  as  t-*0  and  as 
t-»«*  can  now  be  constructed  from  expression  (3.19)  for  the 
Laplace  transform  of  the  ERR.  For  the  short  time 
approximation,  it  is  necessary  to  determine  the  asymptotic 
behavior  of  G(s)  as  s-»<®.  To  this  end,  it  is  required  to 
study  the  behavior  of  mfrs+ai/t)  and  /3=-aeq„  as  s-»».  From 
the  definition  of  the  Carson  transform  m(s) ,  it  is  easily 
seen  that 

2 

Furthermore,  since  iq*  is  the  unique  root  of  G,  (z)=0  in  the 

lower  half-plane  lm(z)<0,  p=p{s,i)  satisfies  the  equation 

~  2  2 
m  ( r  s+oer/J )  =i  ( l  +  (  r  s/aip  )  )  . 


In  the  appendix,  It  is  shown  from  (4.2)  that 

/?--  — 7^ — ' — .  +  o  ( s )  as  s-*«. 

a  [  (c/c,)-* r  ] 

Therefore  if  G(s)  is  written  as  an  asymptotic  series  in 
powers  of  i,  it  is  found  that  as  s-*<* 

S 

G(s)=Jl — -  *-U-T- !  -*  f-  -  ifaltc/c^J-i  ]+®14£Jli,  +  o(s“2)]. 

'  2  TjjT  ( l+€  )  c  [s  r  [  J  2m(  0 )  2  '  _ 


Hence,  if  one  assumes  that  G(t)  has  a  Maclaurin  expansion  in 
a  neighborhood  of  t=0  then  from  (4.4)  and  standard 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 
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At  this  point,  it  is  worth  noting  that  by  allowing  t-,0 
in  (3.19),  one  obtains  the  ERR  for  the  singular  stress 
field,  i.e.  the  energy  flux  into  the  crack-tip  per  unit 
crack  advance  in  the  absence  of  a  failure  zone.  (If  one 
specializes  further  to  steady-state  conditions,  one  obtains 
Atkinson's  expression  for  the  ERR  based  on  a  local  work 
argument  at  the  crack-tip. )  One  then  finds  that  in  this 


(4.6) 


(4.7) 
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r  s 


^=q[lly  j  +  o(s)  for  0<v<c, 


-o°  -(-1/2) 


/,*^['I0rn',<r)  dr] 


i1^2  +  o(s1>/2)  for  v=c*. 


p-po  +  o(l)  for  c*<v<c, 

where  p0  Is  defined  implicitly  through  the  equation 
2 


'o 

2  -« 


-  —  =[  e-a7^°rm' (r)  dr. 

.C#J  Jq 


To  avoid  seperate  cases  in  displaying  subsequent  formulas  it 
is  convenient  to  define  po-0  for  0<v<c„[. 

Again,  from  standard  asymptotic  arguments  for  the 
Laplace  transform 


lim  G ( t ) — 
t-»°° 


Le  ae  ( 1-e 


|  l-(fe/30  >  2  |  (  1/2  > 

TT+T)  ( 1+/J0) 


rn(ai /t  ) 


1  - 


m(ai /t ) 


- ( 1/2 ) 


which  is  precisely  the  steady-state  solution  found  in 


Walton( 1987a) . 

Several  comments  on  these  asymptotic  results  are  in 
order.  First,  it  should  be  observed  that  whether  or  not  the 
model  assumes  a  failure  zone  dramatically  affects  both 
qualitatively  and  quantitatively  the  behavior  of  the  ERR  as 
a  function  of  time,  crack  speed,  and  material  properties. 

For  example,  from  (3.19)  it  is  easily  shown  that  for  any  t, 

L  2a 

if  e>0,  lim  G(t)*»  whereas  if  e=0,  lim  G(t)=— * — -D(t),  where 
7-*0  7_*0  ^ ° 

sD ( 8 } “T+/5TFTS )  and  from  {4,2)  P(s,0)«~[m(rs)  ] (-1/2)  . 

However,  in  the  steady-state  limit  (4.10),  G  approaches  a 

finite  limit  as  the  crack  speed  vanishes,  both  with  and 


(4.8) 


(4.9) 


(4.10) 
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without  a  failure  zone.  Specifically,  it  is  easily  seen 


that  for  the  steady-state  limit 
L  2a 

lim  G(o»)=— g— - £  *2 i  i  for  e>0 
T-»0  H<*> 


Le2S  1 


for  t  =0 . 


Thus  G  becomes  infinite  as  y-*0  except  for  t=0  (no  failure 
zone)  or  under  steady-state  conditions. 

The  reason  for  this  behavior  is  found  through 
consideration  of  the  crack  face  particle  velocity 

«  d  W  A  w 

u, (x, .O.tJs^xtx.O.tJ-v^—  (x,0,t) .  A  consequence  of  the 
o  r  ox 

assumption  that  there  is  an  initial  jump  discontinuity  in 

d  w 

the  applied  crack  face  tractions  is  that  ^(x,0,t)  does  not 

vanish  as  v-*0.  Thus,  from  (3.3)  it  follows  that  lim  G(t)=». 

v->0 

In  contrast,  when  e=0,  one  sees  from  (4.6)  that  G  is  merely 
a  product  of  K(t)*K'(t),  which  depends  only  on  the  SIF  and  a 
simple  function  of  crack  speed  and  glassy  material 
properties  that  is  independent  of  the  crack  face  particle 
velocity  and  that  remains  bounded  as  v-*0.  Moreover,  in 

r°  -  aw 

steady-state,  G  is  given  by  G»-  (x)y^(x,0)  dx  and  thus 

“af 

remains  bounded  as  v-»0. 

Other  differences  between  the  e«0  and  &>0  cases  are 
evident  in  the  short  time  behavior  of  G(t).  In  particular, 
from  (4.5)  it  is  easily  seen  that  when  t>0, 


(4.11) 
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\il  t  ?>o 

whereas,  for  e=0,  it  follows  from  (4.7)  that 

lim  G(t)=0. 
t-*0 

Additionally  from  (4.5),  when  fc>0,  G'(0)  is  given  by 

[a|l-(v/c)]«^l{J|(c./c)]i. 

From  the  fact  that  m'(0)<0  and  using  the  identities 
m' (0)  u ' (0) 

^  q  j  and  a=c„,r/ae,  it  follows  that  the  sign  of  G'(0) 
depends  upon  the  crack  speed  and  material  properties  through 
the  relation 

G ' ( 0 ) >0  Jg'(0)<0  if  and  only  if 

i^|g|l(ae/c)>[l-(v/c)]  [I|l|°|l(ae/c)<(l-(v/c)]]  . 

In  particular  for  fast  enough  crack  speeds,  G'(0)  will  be 
positive,  i.e.  G(t)  will  initially  increase  with  time. 
However,  for  any  given  crack  speed  if  the  combination 
-M^j|^y-L(ae/c)  is  small  enough,  then  G(t)  will  initially 
decrease  with  time. 

Attention  will  now  be  directed  toward  describing  the 
manner  in  which  G(t)  converges  to  the  steady  state  limit  as 
t  approaches  infinity.  It  is  useful  for  comparison  purposes 
to  consider  the  question  first  for  elastic  material.  An 
expression  for  sG(s)  for  elastic  material  emerges  by  letting 
and  m(s)sl  in  (3.19).  Recalling  the  identity 
a*c,r/ae  and  noting  that  c=c„ ,  one  sees  easily  that 


(4.12) 


(4.13) 


(4.14) 
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„f(r,.L.a%  <1-0  L-U/»2l(1/2>  >  , 

w>SiPr  tt«)  (U/.1  1 

where  r t «a  /c.  Moreover,  from  (3.20)  it  is  straight  forward 

© 

to  deduce  that  for  elastic  material 


(4.15) 


/J*sr  ,  /  ( 1-7  )  . 

Substitution  of  (4.16)  into  (4.15)  and  routine  algebraic 
manipulation  yields 
L  2  a 

e  e  (1-4)  ,,  .2J-1/2) 
sG(s)— TRr|  (l-i  ) 

(l»^|l  11/21  I  Hsr.t/i  I  (lt|Ii)"1(l+£lii.)  l‘1/2>  . 

It  is  instructive  to  consider  first  the  elastic  limit 

with  no  failure  zone.  Letting  t-»0  in  (4.17),  one  then  has 


(4.16) 


L  2a 

sG(s)»-|^  e  { 1-7  2  )  * _1/2  )  (  i+^Ll) -1  from  which  it  follows  that 
L  2a 

G(t)=  |  e  (1-72) (_1/2) (i-e'ait)  where  at = ( 1-t  ) /r  t = ( c-v) /a0. 
Thus,  G ( t )  converges  to  its  steady  state  limit  with  an 
exponential  decay  rate  that  decreases  monotonically  with 
increasing  crack  speed  up  to  the  shear  wave  speed. 

A  somewhat  more  complicated  expression  for  G(t)  results 
for  OO.  Calculating  G(t)  by  Laplace  inversion  of  (4.17) 
may  be  conveniently  done  by  defining  a,=(l-7 )/r,  (as  above), 
aj*7/(r,e),  a1=a1/e  and  a^=(l+7 )/(r,e ) ,  and  then  considering 
the  Laplace  inversion  of 

H(s)-i  (l+s/a, ) (1/2) (1+s/a4) (_1/2)  .  It  is 

convenient  to  factor  H(s)  as  the  product  H(s)=H4 (s)H, (s) 

and 


(4.17) 
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H,(s)=7i7a^[[s+(a,+a4)/2]2-[(a,-a4)/2]2]<":1/2).  Then  H(s) 
can  be  inverted  to  h(t)*hj*ha  where  h4 ( t ) =bt a ( t ) +l+ba e  a‘* 

and  I0(li  )  where  b.-Ul*. 

b**-l+(4 /i ) ( 1-* ( l-i ) ) ,  6  { t )  is  the  Dirac  point  measure,  and 
I  (t)  is  the  modified  Bessel  function  of  the  first  kind  of 


order  zero  defined  by  I  ( t )  = 


may  be  written  as 


*  (^)2r 

^  (  r  1  )  2 

r=0  '  ' 


Therefore  the  ERR 


2 

n,  +  s-  e  ae  (I-*-  )  , ,  2  4  (-1/2) 

G(t)— 2^—  TT+D(1  7  } 


b1h2(t)  +  J  ha  { t )  dt  +  b2e~ait*h2(t)] 


From  the  fact  that  0<al<a,<a.,  an  easy  argument  shows  that 

4 

^h(t)  decays  like  e-ait  as  t-*<».  Thus  G(t)  converges  to  its 
steady-state  limit  with  the  same  exponential  rate  as  in  the 
t=0  case. 

It  is  easily  seen  that  for  a  viscoelastic  material  the 
situation  is  considerably  more  complicated  due  to  the 
combined  Influence  of  material  Inertia  and  viscoelastic 
stress  relaxation  upon  convergence  to  steady  state.  A 
general  property  of  the  Laplace  transform  ?(s)  of  a  function 
f(t)  is  that  f(t)  decays  exponentially  in  time,  say 
f(t)~e_at  with  a>0,  if  and  only  if  F(s)  is  analytic  in  the 
halfplane  Re(s)>-a.  The  expression  (3.19)  for  sG(s)  is 
valid  for  s  real.  Its  analytic  extension  for  complex  s  is 
given  by 


(4.18) 


r 
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L  2a 

sG  ( 8 )  -  |  |  (  l+/> )  *  I®(  rs+ar  /t  )  ] 


[!-(*/»)  ]/ 


[— ^—  2  -  ] 

L(mrs/a)  B(r8+ai/t)J 


(1/2) 


Determining  the  largest  value  of  a,  a  ,  for  which 

max 

sG(s)  is  analytic  in  the  halfplane  Re(s)>-a  is  a  difficult 

task  that  clearly  depends  upon  the  particular  details  of  the 

transform  m(rs+a7/t).  From  the  previously  stated  properties 

assumed  for  m(t)  it  follows  that  [m( r s+ai /«. ) ] -1  is  analytic 

for  Re(s)>-^-.  If  m(t)  is  a  powerlaw  in  t,  then  mfrs+air/e) 

has  s=-^l.  as  a  branch  point  and  a  can  be  no  larger  than 

that.  As  a  second  example,  for  a  standard  linear  solid  with 
-t  ~  -1 

m(t)*l+/7e  ,  [m(rs+ai/t)]  is  analytic  for 

Re(s) >-[ {rj  +  1 ) _1+a t /e)/r . 

The  heart  of  the  matter  lies  in  determining  the 
analyticity  properties  of  p( s,y),  which  is  defined 
implicitly  through  equation  (4.2).  Whether  or  not  p{s,y)  is 
analytic  in  some  half plane  Re(s)>-a,  a>0,  depends  upon  the 
particular  way  m(t)  decays  to  its  equilibrium  value  as  t-»«*. 
This  necessitates  a  case  by  case  analysis  for  different 
forms  of  m(t).  We  content  ourselves  here  with  illustrating 
the  differences  that  exist  between  materials  with 
exponentially  decaying  modulus,  such  as  a  standard  linear 
solid,  and  those  for  which  m(t)  decays  as  a  power  of  t  to 
its  equilibrium  value,  such  as  a  simple  powerlaw  material 


with  m(  t  )®l+/7  ( 1+t ) 


-n 


n>0. 


(4.19) 
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An  important  observation  to  be  made  for  powerlaw 

material  is  that  G(t)  cannot  have  exponential  decay  to  its 

steady-state  limit  when  7<l,  i.e.  v<c,  .  The  reason  for  this 

is  that  is  not  analytic  at  s*0  as  can  be  seen  from 

the  following  argument.  It  was  remarked  earlier  that 

p(0,i)~0  whenever  i<l,  in  particular  lim  fi  ( s  ,1  )  /s  — r- 

s-0+  a 

Moreover,  for  powerlaw  material  m(s)  has  s=0  as  a  branch 
point.  If  /J(s,7)  were  analytic  at  s=0  then  the  right  hand 
side  of  equation  (4.2)  would  also  be  analytic  there. 

However,  the  left  hand  side  of  (4.2)  has  s=0  ai  a  branch 
point.  This  contradiction  proves  the  claim. 

On  the  other  hand,  for  l<7<m(0) ^ (i.e.  c,<v<c) , 
since  p(0,7)>0,  equation  (4.2)  defines  p(s,7)  as  an  analytic 
function  in  a  neighborhood  of  s=0.  In  particular,  p(s,7)  is 
analytic  in  a  halfplane  Re(s)>-a,  a>0.  However,  finding  the 
largest  such  a  is  difficult.  For  elastic  material,  the  rate 
of  exponential  convergence  of  G(t)  to  its  steady-state  limit 
corresponds  to  the  negative  real  value  of  s  for  which 
p( s,7)=-i.  For  viscoelastic  materials,  again  a  case  by  case 
analysis  will  be  required  to  determine  precisely  where  the 
singularity  of  sG(s)  with  largest  real  part  will  occur. 

The  remaining  case  i*l  (i.e.  v»c,)  is  easily  handled. 
Indeed,  for  a  general  material,  not  just  powerlaw  or  a 
standard  linear  solid,  equation  (4.2)  admits  no  solutions 
for  s<0.  Thus  p(s,l)  is  not  analytic  in  any  halfplane 
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Re(s)>-a,  a>0  and  G(t)  cannot  converge  exponentially  to  its 
steady-state  limit.  Furthermore,  it  is  not  difficult  to  see 
from  equation  (4.2)  that  the  left  endpoint  of  the  largest 
s-interval  containing  zero  for  which  an  admissible  solution 
/j(s,7)  exists  converges  to  zero  as  *r  approaches  1  from 
above . 

These  observations  made  for  powerlaw  material  whenever 
1< i <m( 0 ) ^ 1 *  (c,<v<c)  are  equally  valid  for  a  standard 
linear  solid.  A  departure  in  behavior  occurs  for  0<7<i.  If 
m(  t )  ^l+rje~t ,  then  m(s) -1* -S  and  (4.2)  can  be  rewritten 

as 

l  +  (l+/7)Z_2f  Z  l2 

— i+z'  '  [z^sj  ' 

where  z*r s+cnr/J ( s , i  )  .  When  s>0,  a  root  z  must  be  sought  that 

is  greater  than  re,  whereas  z<rs  is  required  when  s<0  is 

suitably  near  zero.  An  examination  of  the  graphs  of  the 

functions  on  either  side  of  equation  (4.20)  quickly  reveals 

that  admissible  solutions  exist  for  any  s>0  and  s<0  suitably 

near  zero.  Moreover,  no  solution  exists  for  7  =  1  and  the 

left  hand  endpoint  of  the  largest  s  interval  containing  zero 

on  which  a  solution  exists  tends  to  zero  as  s-»0~.  One 

concludes  from  this  that  for  a  standard  linear  solid,  G(t) 

converges  exponentially  to  its  steady-state  limit  for  0<7<l 

and  l<7<m(0)  ^ 1//2^  but  not  for  7-1.  Moreover,  the  rate  of 

exponential  convergence,  i.e.  a  ,  tends  to  zero  as  7-»i±. 

max 


(4.20) 


must  vanish  as 


It  is  also  easy  to  see  from  (4.2)  that  amax 


-r-»m(  0 ) 


(1/2) 
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5.  Conclusions 

The  principle  contributions  of  this  paper  to  the  study 
of  transient  mode  III  crack  propagation  are  the  solutions 
for  the  displacement  and  stresses  for  general  loadings  and 
general  shear  moduli  and  the  inclusion  of  a  failure  zone 
into  the  model  and  the  calculation  of  the  energy  flux,  G(t), 
into  the  failure  zone.  It  was  then  observed  that 
significant  qualitative  and  quantitative  differences  exist 
in  the  behavior  of  G(t)  as  a  function  of  time,  crack  speed, 
and  material  properties  between  a  model  incorporating  a 
failure  zone  and  one  which  does  not. 

The  question  of  the  rate  of  convergence  of  G(t)  to  its 
steady-state  limit  as  t-*<»  was  also  investigated.  It  was 
observed  that  this  rate  of  convergence  depends  in  a 
complicated  way  upon  the  rate  of  stress  relaxation  and  crack 
speed.  In  particular,  for  a  standard  linear  solid  in  which 
stresses  relax  exponentially  fast,  G(t)  converges 
exponentially  fast  for  all  crack  speeds  except  the 
equilibrium  shear  wave  speed.  Thus  for  crack  speeds  near 
the  equilibrium  shear  wave  speed,  it  is  expected  that 
steady-state  conditions  would  set  in  more  slowly  than  for 
crack  speeds  above  or  below  it.  Also  the  exponential  rate 
of  convergence  is  lost  at  the  glassy  shear  wave  speed.  In 
contrast,  for  power  law  material,  an  exponential  rate  of 
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convergence  of  G(t)  does  not  occur  for  any  crack  speeds  less 
than  or  equal  to  the  equilibrium  shear  wave  speed  whereas 
for  speeds  between  the  equilibrium  and  glassy  shear  wave 
speeds  G(t)  does  converge  to  steady  state  at  an  exponential 
rate . 
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Appendix 

The  behavior  of  p  as  s-*»  is  determined  by  equation 
(4.2).  If  one  lets  h=p/s>0  then  (4.2)  becomes 
m(s[r  4<xirh]  )=i  2  ( l+(r  /a-»h)  )  2  . 


The  left-hand  side  m(s[r+cnrh]  )=m(0)  + 

t  2  „ 


J  Cl 


-(r-knh)sr  ,  ,  .  . 

'  m ' ( r)  dr 


has  the  limit  m(0)  = 


>i  as  s-*<*>.  Therefore  (A.l)  will  be 


satisfied  only  if  h=h  +o(l)  where  h  >0.  As  s-*<»,  (A.l) 

00  00 


becomes 
2 


=i2(l+(r/aih  )  )2 


It  is  easily  seen  that  h  =— n — 7 - — t ,  and  thus 

P=a[(c/f;T-->]  +  o(s)  as  s-~- 

To  find  G ( s )  in  terms  of  powers  of  s  1 ,  one  needs  to 
determine  the  expansions  of  each  of  the  individual  factors. 
First,  one  needs  to  find 


a(rs+a?/fe )=m(0)+ 


J  •'<' 

J  n 


( r s+a7 /€ ) rm , ( r )  dr=m(0)+m4s  1+o(s  1). 


-1 


The  s  term  is  found  by  taking  the  limit  as  s-*«>  of 
00 

s[  e  (r8+a7  /*■  )rmi  (r)  dr.  After  integrating  by  parts, 

J  0 

finds  that  m,*^1  (0)  .  Therefore 
■(rs+ai/t )“m(0)+im' (0)s  1+o(s  1)  and 
15(r.«./.,]-1=5^[l  -  *  o(s-»). 

The  expression  lie'll.1 ’ 


one 


(A.l) 


(A. 2) 


(A. 3) 


can  be  expanded  as 


38 


[l  -  atc/c.-i]^ 


The  final  term 


/  "I  t 

+  o  ( s  )  as  s-*» . 
(a/O2  1 


(r s+a-r/fe  ) 


CM*-2)--1 

m 

(1/2) 


=m(  0) 


T57 
1  + 


_  m'  (0)  1 
m( 0 )  r s 
m '  (0 )  1 


Sm(0)rs 


m(rs+ai /e ) 
-(1/2 


+o(s_1) 


-( 1/2 ) 


+o(s  1 ) . 


If  (A.3)-(A.5)  are  substituted  into  (3.19),  one  obtains 
(4.4)  . 

To  determine  the  behavior  of  p  as  s-*0 ,  one  must  again 
consider  equation  (4.2).  Note  that  the  left-hand  side  of 


(4.2) 

00 

m(  r  s+ai  p  )  =m  (  0 )  +  [  e  ( T  s+ai  p  )  rm  ,  (  r  )  <}r>  1  . 

J  0 

For  0<v<c„,  the  limit  as  s-»0  of  the  right-hand  side  of 

(4.2) ,  i 2 ( 1+rs/aip ) 2 ,  will  satisfy  (A. 6)  only  if  p=PiS+o(s) 

as  s-»0 .  Therefore  as  s-*0,  (4.2)  becomes 

l=-»2(l+r/a-»/jl  )2 

The  solution  of  (A.  7)  for  px  is  easily  seen  to  be  p ,  =g  ^  ^ 

For  c*<v<c,  the  right-hand  side  of  (4.2)  satisfies 

i2(l+rs/a7/J)2>72>l. 

If  p-*0  as  s-»0  then  it  can  be  seen  from  (A.  6)  that  the  limit 
of  m( r s-kxi/3  )  =1  as  s-*0  which  contradicts  (A.  8).  Thus 
P*P0+o(l),  po>0  as  s-*0  where  po  satisfies 

72-(c/c,)2*f  e~ai^°rm' (r)  dr. 

J  0 

For  v=c*.  it  can  be  shown  that  to  satisfy  (4.2)  that 
P-+0  and  s/p-tO  as  s-*0.  Therefore  consider  p-p^  sr}+o(srf), 


(A. 4) 


(A. 5) 


(A. 6) 


(A. 7) 


(A. 8) 
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0 <>7 <  1  ,  as  s-*0.  If  one  substitutes  this  into  (4.2)  and 
simplifies  the  equation,  one  obtains 

-oo  1  -q 

-ay  p ,  s^  J  rm’(r)  dr  +  o(sn)=^~ —  +  ofs1"'7)  as  s-*0 

00 

assuming  that  rm'(r)  dr  exists.  If  this  integral  does  not 
J  n 


exist  then  one  must  do  the  asymptotics  for  the  particular 

shear  modulus  needed.  Since  the  coefficient  of  each  side  is 

(-1/2) 

I’m  '  ( r )  dr  I 

0 


non 


JTr  r  r00  1  (  ~  1  /  2  ) 

-zero,  q  =  l/2  and  /J,= -  -  rm'(r)  dr  .  In  summary, 

’No 


rs 


P  ~a  [  fT'-y  ']  +  °<S)  f°r  0<V<C*< 


—  -  rm ' ( r )  dr 

1  L  J  0 


(-1/2) 


s1/2  +  o(s1/2 )  for  v=c, , 


(A. 9) 


p—p0+o ( 1 )  for  c,<v<c. 


and 
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